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ABSTRACT 



There is currently a debate over the existence of claimed statistical anomalies in the Cosmic Microwave Background 
(CMB), recently confirmed in Planck data. Recent work has focussed on methods for measuring statistical significance, 
on masks and on secondary anisotropics as potential causes of the anomalies. We investigate simultaneously the method 
for accounting for masked regions and the foreground integrated Sachs- Wolfe (ISW) signal. We search for trends in 
different years of WMAP CMB data with different mask treatments. We reconstruct the ISW field due to 2 Micron 
All-Sky Survey (2MASS) and NRAO VLA Sky Survey (NVSS) up to i = 5, and we focus on the Axis of Evil (AoE) 
statistic and even/odd mirror parity, both of which search for preferred axes in the Universe. We find that removing the 
ISW reduces the significance of these anomalies in WMAP data, though this does not exclude the possibility of exotic 
physics. In the spirit of reproducible research, all reconstructed maps and codes will be made available for download at 
http : //www. cosmostat . org/anomaliesCMB .html 



1. Introduction 

In recent years, several violations of statistical isotropy have 
been reported on the largest scales of the Cosmic Microwave 
Background (CMB). A low quad r upole was reporte d in 

and 



COBE data flHinshaw et "aD|l996 - 
confirmed later with WMAP dat 



Bond et al. 
(Sperge 



..199 8 
et al 



2003 



Gruppuso et al.|[2013 ). The octopole presented an unusua 



planarity and its phase seemed correlated with that of the 



quadrupole (Tegmark et al. 


2003 de Oliveira-Costa et al. 


2004 


Slosar STSeljak 2004 C 


opi et al.|2010f. Other anoma- 



2004 | |Bernu i et al 
let al.[ |2004 



C 



ruz 



. 2006 
et al. 12005 



an anoma lous cold spot 

a lignm e nts of other 



2006) 



Vi 



large- scale multipoles QSchwarz et al. 2004| |Copi et al. 
2006), the so-called 'Axis of Evil' (AoE) I pmd fe Magueijo 



2005a ) and other violations of statistical i sotropy (Hajian & 
Souradeep 2003; Land & Magueijo 2005b). Recently, Planck 
has confirmed that these large-scale anomalies are stil l 
present in the CMB (Planck Collaboration et al. |2013a|b I, 
ruling out any origin due to a systematic in the data. 

These anomalies are interesting because they point to- 
wards a possible violation of the standard model of cosmol- 
ogy, which predicts statistical isotropy and Gaussian fluctu- 
ations in the CMB, and ther efore offer a window into exotic 
early-universe physics (e. g., Liu et al. |2013a|b Wang|20"T3 
Buckley & Schlegel 20131. However, there is still much de- 
bate over the possible causes of these anomalies. Since these 
effects are on very large scales where there is large cosmic 
variance, the statistics use d to measure the significance of 
the anomalies are subtle (Bennett et al. 2011 Efstathiou 



et al. 


2010 


Gold et al. 


2011 



be either contamination due to Galac tic foreground effects 
( Hinshaw et al.|20T2 Kim et al.|2012 1 or cosmological fore- 



grounds that lead t o secondary ani sotropies in the CMB 
([Rassat et al.||2007[ |Rudnick et aT1|2007[ [Peiris k, Smith 



120101 ISmith fc Huterer||2TO 



Yershov et al.||2012| |FrancTs 



2013) 



|fc Pe acock 2010; Rass at et al. 

Rassat et alj |2013|) have simultaneously investigated 
asks 



the impact of masks and the integrated Sachs- Wolfe (ISW) 
effect. Missing data were accounted for with the sparse in- 
paintin g technique described in |Starck, Fadili, k, Rassat] 
(2013b), which does not assume the underlying field is ei- 
ther Gaussian or isotropic, but allows for it to be. It was 
found that removing the ISW reduced the significance of 
two anomalies in WMAP data (the low quadrupole and 
the quadrupole / octopole alignment) . In this work we focus 
on two anomalies, both related to preferred axes on the 
sky: the AoE effect and mirror parity. In Section [2] we de- 
scribe the reconstruction of the ISW field from 2MASS and 
NVSS data. In Section [3] we search for violations of statis- 
tical isotropy before and after inpainting, as well as after 
both inpainting and ISW subtraction. In Section |4] we dis- 
cuss our re sults and summaris e our results combined with 
those from Rassat et al. (20131 in Table [T] 



anais. rassat Qepfl.ch 



2. Estimating the large-scale primordial CMB 

2.1. Theory 

Since statistical isotropy is predicted for the early Universe, 
analyses should focus on the primordial CMB, i.e. one 
from which secondary low-redshift cosmological signals 
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have been removed. In Appendix [D] we briefly review how 
to estimate the primordial CMB with a reconstructed map 
of the integrated Sachs- Wolfe effect. In practice, we esti- 
mate the primordial CMB on large scales by: 



<5prim — C>OBS 



?2MASS 
°ISW 



cNVSS 
°ISW 



(1) 



where S 2 g^ ss and Sf^ s are the estimated ISW contri- 
butions from the 2MASS and NVSS surveys (see Section 
[Z2| Appendix |D] and Section 2.2 from |Rassat et al.|[20T3| 
for details on how these are estimated). The terms Sobs 
and S pi i m correspond to the observed and primordial CMB 
respectively. The term SkD,i=2 is the t emperature signal 



due to the kinetic Dop pler quadrupole (Copi et al. [2006 
, for which we have produced pub 



Francis fc Peacock|2010 



licly available maps in Rassat et al. (20131. 



2.2. Data 



We are interested in identifying trends and therefore con- 
sider WMAP renditions from various years, both before 
a nd after inpainting. These maps are described in Table 1 
of|Rassat et al.| ( |2013[ ) (se e also Appendix [Al The 2MASS 
( Jarrett|2004f anT!WSS ( [Condon et al |1998|) data s ets ar e 



described in Sections 3.2 and 3.3 or in Rassat et al. (20131. 



We use the reconstructed ISW effect due to 2 MASS and 

pOlI) and 



NVSS galaxies, usin g the method in Dupe et al 



Rassat et al. ( |2013 l, but reconstructing the ISW fields up 
to £ = 5. The ISW reconstruction method is cosmology in- 
dependent, and estimates the ISW amplitude directly from 
the CMB and galaxy field cross-correlations. This means 
that there is a different ISW map for each CMB map ren- 
dition considered. All map reconstructions are done using 
the publicly available I SAP co de, and the specific options 
used are described in detail in Rassat et al. (2013). 



The large-scale ISW temperature field (I = 2 — 5) due 
to 2MASS and NVSS (where the amplitude is estimated 
from a correlation with WMAP9 data) is plotted in Figure 
[I] Since there is little redshift overl ap between the surveys 
(se 



(see Figure 3 in [Rassat et 



al. 2013), we estimate the ISW 



contribution from each survey independently and add the 
two resulting ISW maps to produce the map in Figure [T| 

In Figure pi we plot the amplitude of the ISW tempera- 
ture power C^ v (£) for £ = 2 — 5, which is measured directly 
from the map plotted in Figure [T] showing the theoreti- 
cal prediction for a 'Vanilla-model' cosmology and taking 
contributions into account from both 2MASS and NVSS 
galaxies. Data points are shown for WMAP9, 2MASS, and 
NVSS data with Gaussian errors bars estimated analyti- 
cally and assuming 

/NVSS = o 66 ( Le _ the sma ii est f sky 

value amongst the three maps). In general, we find that 
the ISW signal is below what is expected from theory for 
the fiducial cosmology we assumed. 



3. Preferred axes in WMAP data before and after 
inpainting and ISW subtraction 

We focus on both the AoE statistic and the even/odd mirror 
parity anomalies, which are described in detail in Appendix 
[B[ In Appendix [C] we validate that sparse inpainting pro- 
vides a bias-free reconstruction with which to study these 
anomalies using a large set of simulated CMB maps. 



Total ISW (2MASS + NVSS) 




-20.0 



20.0 



Fig. 1. The large-scale ISW temperature field (£ = 2 — 5) 
due to 2MASS and NVSS galaxies. The amplitude is es- 
timated directly from the data, including cross-correlation 
with WMAP9 data. 



400 



300 



ST 200 

4 



O 100 

It 



-100 



-200 



Theory (2MASS+NVSS) 
Estimated 




Fig. 2. Amplitude of the ISW temperature power (£) 
(in /J.K 2 ) for I = 2 — 5. The solid line is the theoretical pre- 
diction, and the data points are estimated from 2MASS and 
NVSS galaxies with WMAP9 data are shown with Gaussian 
error bars assuming / s k y = 0.66, i.e. the sky coverage of the 
NVSS survey. 



3.1. Axis of Evil (AoE) 



Rassat et al. (2007) has searched for an AoE directly in 
2MASS data, hoping that this could link the measured 
anomaly to the ISW effect, but has found no preferred axes 
in 2MASS data. Here, we can test an estimate of the pri- 
mordial CMB estimated directly. 

In Table |D.1| we report the preferred modes and axes 
for the 11 CMB maps considered in this paper. This table 
can be directly compared with Table 1 in |Land fc Magueijo| 
(2007) which reviewed the preferred modes and axes for 
Wl and W3 data. By analysing the 'raw' maps (section 1 
of Table D.l I, we confirm their result that the preferred axis 



that was present in first-year data is not present in the ILC 
W3 data. We find the mean interangle is anomalously low 
for 4 WMAP renditions (TOH1, W5, W7, W9) wit h 9 ~ 
20°, but tha t for W 3 data, 8 ~ 52°, similar to what Land 
|fc Magueijo| d2007l found for Wl and W3 maps in their 
updated analysis of the AoE. As they found, the change 
occurs because the preferred mode for £ = 2 is m = 2 
for Wl (and also years 5 to 9), whereas for W3 data, the 
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preferred mode for £ — 2 is m = 0. They attribute this 
fact to the discontinuous nature of the AoE statistic, and 
underline how this feature constitutes a weakness in the 
AoE statistic. 



However, following Francis & Peacock (2010), we sub- 
tract the kinetic Doppler (kDj from the public WMAP 
maps (i.e. by subtracting the last term in Equa tion IT} and 
re-perform the AoE analysis (section 2 of Table D.l| where 
only the quadrupole 1 = 2 has changed) . We find the statis- 
tic is now stable across all WMAP data considered, with 
the mean interangle spanning 8 = 18.5 — 20.8°. 

The results after inpainting of the CMB maps (and kD 
subtraction) are presented in section 3 of Table D.l| We 



find inpainting has no effect on the preferred modes and di- 
rections for £ = 2 and £ = 5, which are generally unchanged 
and are the same for all maps, as was the case before in- 
painting. For £ = 4, the preferred mode and direction for 
all WMAP renditions now becomes similar to that for the 
TOH map before inpainting. The octopole (I = 3) is the 
scale that changes the most after inpainting, both preferred 
mode and direction are changed for all WMAP renditions. 
The preferred modes and directions for £ = 2 — 5 are quite 
stable across all WMAP renditions after inpainting. We 
note t hat when searching for quadrupole/octopole align- 



ment, Rassat et al. ( |2013[ ) showed that the preferred axis 
of the octopole was stable after inpainting. However, this 
statistic enforced a search within a planar mode (m = 3), 
which explains the difference in stability for the octopole. 

The change of only one preferred mode, notably the 
octopole in the W9 rendition, that prefers the planar mode 
m = 3 rather than m — 1 for other WMAP renditions, has 
a significant effect on the mean interangle 9, that changes 
from 48.3 - 51.0° for Wl-7 data to 16.3°, a value which 
is only found in 0.1% of simulations. Further studies with 
improved (and smaller) masks, either on WMAP or Planck 
data, could point to why the AoE statistic on W9 data is 
so different. 

The results after inpainting and ISW subtraction are 
presented in section 4 of Table |D.1| For the quadrupole, 
the preferred mode remains the same, and the preferred 
direction changes slightly. For £ = 4, 5 there is virtually 
no change in either the preferred mode or direction, which 
can be explained by the low values of the ISW signal on 
these scales (see Figure 2| . The scale with the most change 
is the octopole, where both preferred mode and direction 
change for all WMAP renditions. In general, there is good 
agreement across different WMAP renditions, except some 
differences for the octopole, e.g. for W3 and W9 renditions. 
After inpainting and ISW subtraction, the mean interangle 
now varies between 33.1 — 47.8, which corresponds to 1.9 — 
13.3% of values found in simulations. Even for the lowest 
value of 1.9%, the original anomalous alignment is no longer 
as significant as the previous value of 0.1% on the ILC maps. 

3.2. Parity 

We calculate the S± values for each of the WMAP ren- 
ditions using Equation |B.3| and results are presented in 
Table |D.2| where the significance is calculated using 1000 
full-sk y Ga ussian realisations. Before inpainting (part 1 of 
Table D.2|, we find no significant even (5+) or odd 
mirror parity, except for the TOH map, which shows a mild 
preference for even mirror parity (only 2.6% of simulations 
have a higher S+ value). After sparse inpainting (part 2 



of Table D.2|, we find an increase in both even and odd 
mirror parity anomalies (except for the TOH map), with 
only 3.2 — 3.6% of simulations having a higher value for 
S- than both W7 and W9 maps and only 0.90 - 3.1% of 
simulations having a higher value for S+. The significance 
of these remains at < 3cr though. We find the mirror-parity 
anomalies do not persist after sparse inpainting and sub- 
traction of the reconstructed ISW signal (part 3 of Table 
D.2|, meaning the ISW could explain these anomalies. 



4. Discussion 

One of the main successes of the standard model of cosmol- 
ogy is its prediction of the Gaussian random fluctuations 
observed in the CMB. However, even since COBE data, sev- 
eral signatures of lack of statistic isotropy, or 'anomalies', 
have been reported on large scales in WMAP data and 
recently confi rmed in Planck data (Planck Collaboration 
et al.|2013a|b I. Recent focus has been on testing the impact 
of different reconstruction methods or methods for dealing 
with Galactic foregrounds, while others have investigated 
how various foreground cosmological signals could affect 

these anomalies. 

In |Rassat et al. (20131, we found 
ing the ISW signal dueto 2MASS 



that subtract- 
and NVSS data 



due to 

from WMAP data lowered the significance of two pre- 
viously reported anomalies: the low quadrupole and the 
quadrupole/octopole alignment. The octopole planarity 
was also investigated and found to be less significant in 
later WMAP years. 



Anomaly 



From Rassat et al. 
Low Quad 

Quad/Oct alignment 
Oct planarity 



After sparse 
inpainting 



After ISW 
subtraction 



More anomalous 
Less anomalous 
Less anomalous 



Not anomalous 
Not anomalous 
Not anomalous 



This work 
Axis of evil 
Even mirror parity 
Odd mirror parity 



Less anomalous 
More anomalous 
More anomalous 



Not anomalous 
Not anomalous 
Not anomalous 



Table 1. Summary of results in this paper and Rassat et al. 
(|2013[ ). The anomalies in WMAP could be explained by the 
ISW effect, though other explanations remain possible. 

In this work, we continued investigation of two other 
anomalies, both related to preferred axes in the sky: the 
Axis of Evil (AoE) and even/odd mirror parity in CMB 
data, i.e. parity with respect to reflections through a plane. 
We first investigated whether sparse inpainting can be con- 
sidered a bias- free reconstruction method for the two statis- 
tics and found that this is the case. We then applied sparse 
inpainting on various CMB maps up to £ = 5 and recon- 
structed the ISW maps from 2MASS and NVSS data also 
up to I = 5. We consider the significance of the two re- 
ported anomalies, before inpainting, after inpainting, and 
both after inpainting and ISW subtraction. 

Our first approach to the AoE was to remov e the kD 
quadrupole (following Francis fc Peaco'ck] [2010), and we 
found that the AoE is consistent across all renditions of 
WMAP data (TOH W3 , W5, W7 and W9), unlike what 
Land & Magueijo (20071 found. After sparse inpainting, 
we found the AoE is no longer anomalous, mainly due to 
the change in preferred mode and axis of the octopole, 



3 



A. Rassat and J.-L. Starck: Preferred Axes in WMAP CMB Data 



except for WMAP9 data, where the anomaly persists. 
Further studies with improved (and smaller) masks, either 
on WMAP or Planck data, could point to why the AoE 
statistic on WMAP9 data is so different. After sparse in- 
paiting, both even and odd mirror parities are increased in 
significance, but not enough to be considered significantly 
anomalous. 

We found that subtraction of the ISW effect due 
to 2MASS and NVSS galaxies, reduces the significance 
of these anomalies. These results, along with those in 
Rassat et al. (2013) relating to the low quadrupole, the 
quadrupole/octopole alignment and the octopole planarity 
are summarised in Table [T] We note, however, that there 
are other signatures of statistical anomalies on large scales 
that we have not tested (e.g.: north/south asymmetry, cold 
spot, etc.) and that exotic physics remain possible. These 
results are based on WMAP data alone, and should be re- 
peated with Planck. 

In the spirit of reproducible research, all reconstructed maps 
and codes that constitute the main results of this paper will 
be made available for download at http : / / www . cosmostat . org/ 
anomaliesCMB . html 
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Appendix A: WMAP Cosmic Microwave 
Background Maps Considered 



As we are interested in identifying trends in the data, we con- 
sider a suite of 11 different renditions of WMAP data: the 



|Tegmark et all ( j2003 1 reduced-foreground CMB Map (TOH1) 
the Internal Linea r Combination (ILC ) WMAP maps fro m the 



3rd year (IL C W3, [Hinshaw ct al. 2007) , 5th year (I LC W5, [Go!d 
etld~|[2009| 7 th year (ILC W7, |Colcl et al.|[MTl| ), and thelJtn" 



4 http://map.gsfc.nasa.gov 

5 http: / /heasarc . gsfc.nasa.gov/W3Browse/all/nvss . 



html 



year (ILC W9, Hinshaw et al.|2012| ), as well as sparsely inpainted 
versions of these maps. We also consider the sparsely inpainted 
WMAP I LC 5th year data reconstructed by | Delabrouille et al.| 



(De la W5, 2009|l using w avelets. These are summarised in Table 
1 oflRassat et all (12013 1. 
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Appendix B: Statistical anomalies & impact of 
sparse inpainting 

B.l. Preferred axis for low multipoles: Axis of Evil 



It was first noted by de Oliveira-Costa et al. (20041 that both 
the quadrupole and octopole of the CMB appeared planar (i.e. 
anomalously dominated by m = ±£ modes) and were also 
align ed along a similar axis. 

|Land fc Magueijo] ( |2005a[ ) suggest searching for a more gen- 
eral axis by considering the power in any mode m, instead of 
focussing on planar modes. This can be quantified by consider- 
ing their statistic: 



re = max 



Ci m (n) 
•a {2£+l)C e 



(B.l) 



The expressions C^ m (n) are given by CVo(n) = |a«)(n)| 2 and 
Ctm(&) = 2ja^ m (n)| 2 for m > and (21 + l)C e = £ m \a lm \ 2 , 
where ae m (ii) corresponds to the value of the ae, n coefficients 
when the map is rotated to have n as the z-axis. The above 
sta tistic finds both a preferred axis n and a preferred mode m. 

|Land fc Magueijo ( |2005a| ) find that the preferred axes for 
£ = 2, 5 for WMAP 1 data seemed aligned along a similar axis 
in the direction of (£, b) ~ (—100°, 60°), where the I varied from 
~ [—90°,— 160°] and b varied from ~ [48°, 62°]. By considering 
the mean interangle 6 (i.e. the mean value of all possible angles 
between two axes and n^/ for £, £' = 2, 5 and £ £'), they 
find that only 0.1% of simulations had a lower average value 
than the one measured in WMAP1 data (~ 20°) and therefore 
rejected statistical isotropy at the 99.9% confidence level. The 
preferred axis has been dubbed the "Axis of Evil" (AoE). 



B.2. Mirror parity 

Another test is mirror parity in CMB data, i.e. parity with re- 
spect to reflections through a plane: x = x — 2(x ■ n)n, where 
n is the normal vector to the plane. Since mirror parity is yet 
another statistic for which preferred axes can be found (i.e. the 
normal to the plane of reflection), it is complementary to the 
search for a preferred axis described in Section |B.1| 

In pr actice , mirror parity and preferred axes found using 
Equa tion B.l are statistically independent (Land & Magueijo 
20051, so the coincidental presence of both increases the signifi- 



cance of the preferred axis anomaly. With all-sky data, one can 
estimate the 5*-map for a given multipole by: 



(B.2) 



Positive (negative) values of Se(n) correspond to even (odd) 
mirror parities in the n direction. The same statistic can also 
be considered summed over all the low multipoles one wishes to 
co nsider (e.g. the multipole s that have similar preferred axes) as 



Ben-David et al.| ( |2012[ ): Stot(n) = £^1* &(n). The P arit y 
estimator is redefined as S(n) = Stot(n) — (^ max — 1), so that 
(S) — 0. The most even and odd mirror-parit y directions for a 



given map can be considered by estimating (Ben-David et al. 
20121 



max(S) - fi{S) 
a(S) 



and S- 



|min(5) -fi(S)\ 
a(S) 



(B.3) 



where fi(S) and cr(S) are the mean and standard deviation. 

Other s have also studied poin t parity with different statis- 
tics, e.g. Land fc Magueijo] (|2005[) who did not find significant 



point parity in the first WMAP data release. |Kim & Na selsky] 



( 2010a|b| |2011 1, however, find evidence for odd point parity in 
later WMAP renditions, and link this anomaly with the low 
level of correlations on the largest scales. 



Appendix C: Validation of sparse inpainting to 
study large-scale anomalies 

One can test for preferred axes directly on different renditions 
of WMAP data, for example, on Internal Linear Combination 
(ILC) maps. However, these may be contaminated on large scales 
due to Galactic foregrounds ( Hinshaw et al.|2012 I . Another ap- 
proach is to use a different basis set than spherica l harmonics, 
i.e. use a ba sis that is orthonormal on a cut sky (e.g. Rossmanith 



et al. 20121. Alternatively, one can u se sparse inpainting tech"^ 



mques to reconst ruct full-sky maps ( |Plaszczynski et al. |2012 



|Dupe et al.|201 T] |Rassat et al.||2013[) or other inp ainting meth- 
ods, such as diffuse inpainting ( j/acchei et al.|2011|) or inpaintin, 
using const rained Gaussian realisations (|Bucher fc Louis|2012 



Kim et al. 2012J). Any inpainting technique should be teste- 
for potential biases, specifically for the masks and statistical 
tests one is interested in. H ere we use the spar se in painting 



techniques first de scribed in Abrial et al.| (|2~008|) and Starck^ 
Murtagh, fc Fad ili (2010) and refined in |Starck, Fadili, fc Rassat 
(2013d| to reconstruct regions of missing data. The advantage 
of this method is that is does not assume the 'true' map i s either 
Gaussian or isotropic, yet it allows it to be (see for e.g. Starck 



et al.|2013a) 



|Rassat 
bias- free 



et al. 



(20131 show sparse 
method for the 



inpainting is a 
low quadrupole, 



reconstruction 

quadrupole/octopole alignment, and octopole planarity tests 
Here, we test whether sparse inpainting is a bias-free method 
for both the AoE and mirror parity. We calculate 1000 Gaussian 
random field realisations of WMAP7 best-fit cosmology using 
the WMAP7 temperature analysis mask (/sky = 0.78). 

C.J. Recovering the mean interangle (6) with a realistic 
Galactic mask 

We com pare the mean interangle 6 for the statistic given in 
Equation |B.l| for £ = 2 — 5, for each map before and after inpaint- 
ing, using nside = 128 for the CMB maps. We find 9 ~ 57° ±9° 
in our simulations before and after inpainting, i.e. what is ex- 
pected in the case of isotro pic axes and a Gaussian random field 
( |Land fc Magueijo|[2005a[ ). After inpainting with the WMAP7 
mask, we find that (#true — #in P ) ~ —0.55° ±10.7°, showing there 
is no significant bias in the estimation of the mean interangle af- 
ter sparse inpainting is applied. While the bias is small, the error 
bar on the mean inte rangle after sparse inpainting is n ot negli- 
gible (10°). Following [Starck, Fadili, fc Rassat] ( |2013b[ ), we test 
the statistic on an optimistic Planck-like mask with f B ^ y = 0.93 
and find (0 tT uc — #m P ) ~ 0.17° ± 7.3°, showing we can expect 
better reconstructions with Planck data and smaller masks. 



True 



After Inp. 



Bias 



57.5° ±9.2° 57.0° ±9.2° 



0.55° ± 10.7° 



S- 



2.59 ±0.30 
2.81 ±0.35 



2.59 ±0.30 
2.82 ±0.35 



-0.00039 ± 0.22 
-0.0049 ±0.26 



Table C.l. Mean interangle 8 for the 'Axis of Evil' statistic 
and even (+) and odd (— ) mirror parity statistics S± be- 
fore and after sparse inpainting on 1000 Gaussian random 
field realisations of CMB data and using the WMAP7 tem- 
perature analysis mask for the inpainted maps. The bias is 
taken by considering the difference (true — inp). 

C.2. Recovering mirror parity statistics (S±) with a realistic 
Galactic mask 

To test for possible biases in the S± statistics after sparse in- 
painting, we calculate S+ and S- for each CMB simulation be- 
fore and after inpainting, setting nside — 8 for the CMB maps, 
and nside = 64 for the parity map s (calculated usin g Equation 
B~2| , as in |Ben-David et al.|p(5T2| . As in Figure 6 of |Ben-David 



5 



A. Rassat and J.-L. Starck: Preferred Axes in WMAP CMB Data 



et al. (20121, we find that the distributions of S+ and S- pop- 
ulations do not cha nge b efore and after inpainting ('True' and 



'After Inp.' in Table C.l I . We do not find any signi ficant bias in 



the S± meas urements ('Bias' column in Table C.l I. 

Following Sta rck, Fadili, fc Rassatj ( |2013b[ ), we also test the 



statistic on an optimistic Planck-like mask with / s k y = 0.93 
and find AS+ = -0.0021 ± 0.081 and AS- = 0.00091 ± 0.10, 
showing we can expect significantly better reconstructions with 
future Planck data and smaller masks. 



Appendix D: Recovering the primordial CMB 

Since statistical isotropy is predicted for the early Universe, 
analyses should focus on the primordial CMB, i.e. one from 
which secondary low-redshift cosmological signals have been re- 
moved. The observed temperature anisotropics in the CMB, 
<5obs, can be described as the sum of several components: 

Sobs = <5 pr im + <$isw on large scales. (D.l) 

where 5 pr im are the primordial temperature anisotropies, and 
<5isw the total secondary temperature anisotropies due to the 
la te-time Integrated Sachs Wolfe (ISW) effect (see Section 2.1 
of |Rassat et~aL|2013 1 . 



In practice, the temperature ISW field can be recon structed 



in spherical harmonics, 5}™ , from the LSS field ae m ( Boughn 
et al.|1998||Cabre et al.|2007||G~iannantonio et al.[2008 iT 



<-ISW _ Cg T {t) m „s 

where ge m represent the spherical harmonic coefficients of a 
galaxy overdensity field g(8, <j>). 

The spectra C gg and C 9 t are the galaxy (g) and CMB (T) 
auto- and cross-correlations measur ed from the data or their 
theoretical values (see Section 2.2 of |Rassat et al.|2013 |. 



Map 



1) Before inpainting 










TOH1 


3.25 


(2.6%) 


3.15 


(16%) 


WMAP3 


2.88 


(17%) 


2.81 


(48%) 


WMAP5 


2.93 


(14%) 


2.85 


(43%) 


WMAP7 


2.93 


(14%) 


2.88 


(39%) 


WMAP9 


3.00 


(9.9%) 


2.93 


(34%) 


2)After inpainting 










TOH1 (inp) 


3.09 


(5.7%) 


3.29 


(10%) 


WMAP3 (inp) 


3.31 


(1.9%) 


3.39 


(6.5%) 


WMAP5 (inp) 


3.30 


(2.0%) 


3.40 


(6.1%) 


W5 Dela (inp) 


3.41 


(0.90%) 


3.57 


(3.1%) 


WMAP7 (inp) 


3.34 


(1.4%) 


3.55 


(3.2%) 


WMAP9 (inp) 


3.20 


(3.1%) 


3.54 


(3.6%) 



3)After inpainting 
and ISW subtraction 

TOH1 (inp) - ISW 2.87 (17%) 3.21 (13%) 

WMAP 3 (inp) - ISW 2.94 (13%) 3.20 (14%) 

WMAP 5 (inp) - ISW 3.05 (7.5%) 3.26 (11%) 

W5 Dela (inp) -ISW 3.01 (9.3%) 3.25 (11%) 

WMAP 7 (inp) - ISW 3.01 (9.3%) 3.37 (7.3%) 

WMAP9 (inp) - ISW 2.92 (14%) 3.40 (6.1%) 

Table D.2. Values of even (5+) and odd (S-) parity scores 
for 2 < I < 5 for WMAP data from different years before 
(1) and after inpainting (2), and after subtraction of the 
ISW effect due to both 2MASS and NVSS galaxies (3). The 
occurrence for 1000 full-sky Gaussian random simulations is 
given in brackets. The kD quadrupole has been subtracted 
for all maps. 
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Map 


Mean 




= 2 






= 3 






= 4 




i- 


= 5 






interangle 9 (°) 


(6,0 


m 


(6,0 


m 


(6,0 


m 


(i 


.,o 


m 


Before inpainting 


1) 

TOHl 


20.9 


(58.9, 


-103.4) 


2 


(61.9, 


-104.8) 


3 


(57.8, 


-164.0) 


2 


(47.7, 


-132.6) 


3 


W3 


51.9 


(-27 


.4,3.3) 





(62.3, 


-103.8) 


3 


(34.6, 


-132.2) 


1 


(47.4, 


-129.9) 


3 


W o 


in 7 


(01.2, 


— Izl. t ) 


z 


(62. ,3, 


— lUo.o ) 


Q 


(34.2, 


— lol.o ) 


± 


//I 7 /I 


— Izy.y J 


Q 
O 


W7 


20.4 


(62.7, 


-123.4) 


2 


(62.7, 


-104.0) 


3 


(33.9, 


-131.5) 


1 


(47.4, 


-130.7) 


3 


W9 


19.1 


(60.1, 


-120.6) 


2 


(62.7, 


-105.2) 


3 


(34.2, 


-131.1) 


1 


(47,7, 


-130.2) 


3 


2) 

TOHl - kD 


20.8 


(56.6, 


-106.5) 


2 




















W3 - kD 


20.3 


(62.7, 


-129.5) 


2 




















W5- kD 


18.9 


(57.8, 


-125.7) 


2 




















W7 - kD 


19.5 


(58.5, 


-127.6) 


2 




















W9 - kD 


18.5 


(57.0, 


-124.9) 


2 





















After inpainting 



3) 



TOHl (inp) - kD 


51.0 


(48.5,-116.4) 


2 


(29.3,82.6) 


1 


(57.8, 


-168.1) 


2 


(47.4, 


-140.0) 


3 


W3 (inp) - kD 


48.3 


(59.7,-140.0) 


2 


(31.7,81.6) 


1 


(58.2, 


-165.9) 


2 


(47.0, 


-135.0) 


3 


W5 (inp) - kD 


49.0 


(54.7,-135.0) 


2 


(31.4,81.9) 


1 


(58.2, 


-165.9) 


2 


(47.0, 


-135.0) 


3 


Dela W5 (inp) - kD 


49.0 


(54.0,-140.6) 


2 


(28.0,82.6) 


1 


(58.9, 


-166.6) 


2 


(46.2, 


-135.8) 


3 


W7 (inp) - kD 


48.5 


(55.5,-138.9) 


2 


(30.7,82.6) 


1 


(58.2, 


-165.9) 


2 


(47.0, 


-135.8) 


3 


W9 (inp) - kD 


16.3 


(54.7,-135.0) 


2 


(57.8,-116.4) 


3 


(58.2, 


-165.9) 


2 


(47.4, 


-135.4) 


3 


4) 

TOHl (inp) - kD - ISW 


33.1 


(79.0,-175.5) 


2 


(26.6,-176.8) 


2 


(57.8, 


-168.1) 


2 


(47.4, 


-137.0) 


3 


W3 (inp) - kD - ISW 


35.2 


(75.7, 148.8) 


2 


(26.3,-177.2) 


2 


(58.2, 


-166.9) 


2 


(47.0, 


-135.0) 


3 


W5 (inp) - kD - ISW 


47.8 


(74.6,161.8) 


2 


(20.1,-76.3) 


2 


(58.2, 


-166.9) 


2 


(47.0, 


-135.0) 


3 


Dela W5 (inp) - kD - ISW 


32.6 


(72.0,162.6) 


2 


(27.0,-177.9) 


2 


(58.5, 


-167.8) 


2 


(46.2, 


-135.8) 


3 


W7 (inp) - kD - ISW 


33.4 


(74.2,158.0) 


2 


(27.0,-177.2) 


2 


(57.8, 


-167.1) 


2 


(47.0, 


-135.8) 


3 


W9 (inp) - kD - ISW 


44.1 


(75.0,161.3) 


2 


(54.7,45.9) 


2 


(57.8, 


-167.1) 


2 


(47.7, 


-135.8) 


3 



Table D.l. Preferred axes for multipoles I = 2 - 5 for different WMAP CMB maps for nside = 128. 



7 



